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ABSTRACT 

Observations show significant Ha-emissions in the Galactic halo near 
the edges of cold gas clouds of the Magellanic Stream. The source for 
the ionization of the cold gas is still a widely open question. In our 
paper we discuss the critical velocity effect as a possible explanation for 
the observed Ha-emission. The critical velocity effect can yield a fast 
ionization of cold gas if this neutral gas passes through a magnetized 
plasma under suitable conditions. We show that for parameters that are 
typical for the Magellanic Stream the critical velocity effect has to be 
considered as a possible ionization source of high relevance. 

Subject headings: Galaxy: Magellanic Clouds - Galaxy: halo - ISM: 
clouds 



1. Introduction 

The Magellanic Stream (MS) consists of a 
narrow band of high- velocity neutral hydro- 
gen gas (HI) which stretches over almost 100" 
on the sky, trailing the Magellanic Clouds in 
their orbit around the Milky Way (Mathew- 
son, Cleary & Murray 1974; Mirabel, Cohen 
& Davies 1979). The Stream was observed 
at 21 cm. Thereby a chain of clouds con- 
nected by a gas of lower density was iden- 
tified. The clouds have been labeled MS I 
through MS VI (Mathewson, Schwarz & Mur- 
ray 1977). They are characterized by a high 
density concentration with a relatively sharp 
gradient on the leading edge, which is defined 
by the direction of the proper motion of the 
Large Magellanic Cloud (Jones, Klemola & 
Lin 1994). The Stream has no known stel- 
lar component (Mathewson et al. 1979) and 
its distance to the Milky Way is unknown, al- 
though its leading end (MS 1) connects to the 
Magellanic Clouds and is presumably located 
at a distance of 50- 60 kpc. Moore and Davis 
(1994) estimate the distance of the tip of the 
Stream (MS VI) to be 20 kpc, whereas Gar- 
diner, Sawa & Fujimoto (1994) give a distance 
value of 60 kpc. 

Bregman (1979) considered the MS as the 
hydrodynamic result of the tidal interaction 
between the Magellanic Clouds and the gaseous 
halo of the Milky Way (Murai & Fujimoto 
1980; Lin & Lyndcn-Bell 1982). This tidal 
scenario is supported by an extension to the 
model (Moore & Davis 1994) which was able 
to reproduce the kinematics of the Stream by 
including the stripping associated with an ex- 
tended ionized disk and drag due to a diffuse 
halo. Nevertheless, the models for the ori- 
gin of the MS remain controversial (Wakker 
& van Woerden 1997; Murah 2000), so that 
the inferences about halo gas properties are 



correspondingly uncertain. 

What do we know from observations about 
the characteristic values for the density and 
temperature of the MS? 

The northern part of the MS (MS III - MS 
MVI) has been surveyed with high sensitivity 
by Mirabel, Cohen & Davies (1979) in order 
to study the small scale structure. Five re- 
gions were investigated, which are character- 
ized by narrow filaments running parallel to 
the main direction of the Stream, and small 
cloud-like concentrations. These clouds have 
typical masses of 15 D'^Mq (here D is their 
distance in kpc). The authors conclude that 
the clouds appear to be unstable over the life- 
time of the Stream as a whole and some exter- 
nal confinement mechanism seems to be nec- 
essary to preserve the observed small scale 
structure. This could be provided by a hot 
plasma of density ~ 4 ■ 10~^ cm~^ at a tem- 
perature of 10^ K. 

Weiner & WiUiams (1996) detected faint 
Ha-emission from several points of the MS, 
i.e. the leading edges of the HI clouds MS II, 
MS HI and MS IV. The emission is strongest 
(emission measures EM: 0.5-1 cm~^ pc) on the 
sharp leading-edge density gradient. They 
suggest that ram pressure from the halo gas 
plays an important role in stripping the Stream 
out of the Magellanic Clouds and they sug- 
gest the existence of a relatively dense {uh ~ 
10~^cm~^) and almost fully ionized gas with 
a temperature of ~ 1.7 x 10^ K in the Galactic 
halo at a radius of 50 kpc. 

Weiner and Williams (1996) rule out pho- 
toionization for the Ha-emission. 

They argue that the HI column densities in 
all observed fields are optically thick to ioniz- 
ing radiation which means that any uniform 
ionizing flux should produce ionization rates 
and Ha emission which are roughly constant 
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from field to field, which is not observed. The 
low number of required UV-photons (10^^ s~^) 
could, in principle, be provided by a few OB 
stars located ~ 100 pc away from the ob- 
served field. However, such stars that would 
be easily detectable by the associated lumi- 
nous Hll-region are not observed. 

Murali (2000) recently discussed some the- 
oretical implications. He reexamined the in- 
teraction of the MS with the surrounding 
medium at large distances in the halo of the 
Galaxy. He was able to show that heating 
dominates over drag forces. Due to their high 
relative velocities (> 200 km s~^) the Hl- 
clouds should evaporate if the density of the 
halo was too high. Thus, he concluded a limit 
for the halo density at 50 kpc to be less than 
lO-^cm-'l 

The strength and orientation of the mag- 
netic field in the halo of our Galaxy is still a 
matter of debate (e.g. Vallee 1998). Accord- 
ing to Beuermann, Kanbach & Berkhuijsen 
(1985) the scale height of the field is ~ 6kpc. 
The distance of the MS is several tens of kpc. 
Given a disk magnetic field of about 6/iG and 
taking into account local field compression by 
the MS clouds the magnetic field in the MS 
is expected to be of the order of several fiG. 
Amplifications of magnetic fields by infalling 
gas clouds is well-known from high velocity 
clouds (e.g. Crutcher 1999; Vallee 1998). 

The interaction of the high velocity clouds 
(HVC) of the MS and the medium of the 
galactic halo is an example for the interaction 
of a fully ionized plasma and a cold weakly 
ionized gas. Such a situation has been con- 
sidered by Alfven (1954) in connection with 
his model for the origin of the solar system. 
He proposed that when a neutral gas streams 
through a plasma across magnetic field lines 
at a velocity that exceeds a certain value Vc a 



discharge like process can occur in which the 
neutral gas component is rapidly ionized. The 
critical velocity is proportional to the ioniza- 
tion potential and inversely proportional to 
the square root of the mass of the neutral 
atom. This effect is usually considered in the 
context of plasmas in which the ion-neutral 
collision rate is negligible. This is exactly 
the case in the interaction regions of the MS 
Hl-clouds and the galactic halo. Given the 
problematic issue of lacking ionization sources 
for the Ha, the critical velocity effect (CVE) 
presents an interesting fundamental plasma 
process that is capable of providing the re- 
quired ionization. The CVE occurs when a 
neutral gas streams through a plasma across 
magnetic field lines with a velocity that ex- 
ceeds a critical value measured by the ion- 
ization potential of the neutrals. The kinetic 
energy of the flow is partly converted to heat- 
ing the electrons which can efficiently ionize 
the neutral gas component. 

It is the aim of our contribution to dis- 
cuss the CVE in the astrophysical system 
MS-Galactic halo as an alternative ionization 
mechanism far away from standard ionization 
sources hke stars. 

In Sect. 2 we review the basic condi- 
tions for the CVE to work. We give lim- 
its for the magnetic field strength, neutral 
gas, plasma density, and plasma temperature. 
Sect. 3 presents the application to the ioniza- 
tion of the MS cloudlets, followed by Sect. 4 
which contains the discussion and some con- 
clusions, in particular, on the general appli- 
cability of the CVE-process in astrophysical 
plasma-neutral gas interactions. 

2. The Critical Velocity Effect 

The interaction of a plasma and a neutral 
gas in the presence of a magnetic field can 
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result in the ionization of the neutral compo- 
nent if the relative velocity component across 
the magnetic field Vrei exceeds a critical value 

Vrel > Vc = e (1) 

V m„ / 

where e, and m„ denote the elementary 
charge, the ionization potential of the neu- 
trals, and the neutral mass. This is the so- 
called critical velocity effect (CVE) (Alfven 
1954; Raadu 1978, 1981; Petelski 1981). The 
efficiency is parameterized by e. Formisano 
et al. (1982) have shown on the basis of 
quasilinear calculations that the efficiency can 
drop to e ^ 0.025 if the ionized neutrals, 
i.e. the newly born ions, have sufficient time 
to isotropize in the form of a ring-like ve- 
locity distribution in the plane perpendicu- 
lar to the magnetic field. The rather poor 
efficiency, or equivalently, the necessary very 
high relative velocities, can be avoided for two 
cases. First, if the ionization frequency is 
much higher than the ion plasma frequency 
(Formisano et al. 1982), and second, if the 
growth rate of the excited lower-hybrid drift 
instability '-fihd that heats the ionizing elec- 
trons is sufficiently high (Brenning 1985), i.e. 
lihd > ( fij = eB/niiC is the ion gyro 
frequency; e, B, rrij, and c denote the elemen- 
tary charge, the magnetic field strength, the 
ion mass, and the speed of light, respectively). 
The first alternative is not of interest in the 
present context and, in fact, should only be of 
importance in dense laboratory plasmas. The 
latter condition can be written as (cf. Bren- 
ning 1985) 

/Tne\i ^^Ue ^ c / me(l + /3e) y 

where Ue = {Annee^ /rrig)^/'^ [rie is the elec- 
tron particle density) and Vie = eB/uieC de- 
note the electron plasma and gyro frequencies 



and (3e = SimekTe/ B^ [k and Tg are the Boltz- 
mann constant and the electron temperature) 
is the electron plasma beta. For high-/? plas- 
mas efficient ionization is possible for electron 
temperatures kT^ > \mivl^i. For low- (3 plas- 
mas the relative flow velocity v^ei has to be 
sub-Alfvenic, Vrei < va = B / {A-KniimiY/'^ , due 
to the second inequality. Here and in the fol- 
lowing we assume quasineutrality rie = rii = 
n. In the low-/? case Eq.(0) can be reformu- 
lated 

1 

7 G-^cm-t < — < 7 — G-^cm-i, (3) 

B Vrel 

if we choose the ion and neutral masses equal 
to the proton mass, i.e. we deal with a par- 
tially ionized hydrogen plasma, or 

Va 

va<.c< c. (4) 

Vrel 

Fig. 1 illustrates the constraints for the mag- 
netic fleld strength and the particle density 
that guarantee efficient heating for Vrei = 
10~^c and for Vrei = 10~^c, for example. In 
the first case, the region between the solid 
line and the dotted one represents the combi- 
nations of particle density and magnetic field 
that allows for the CVE to operate efficiently. 
In the latter case {vrei = lO'^^c), the same 
holds for the region between the solid line and 
the dashed one. 

Up to now we concentrated on the hydro- 
gen component which is dominant in the in- 
terstellar and intergalactic medium. In the 
context of the 'missing ionized helium' puzzle 
(Domgorgen 1997; Reynolds & Tufte 1995), 
it is of interest to discuss the infiuence of the 
CVE on the helium component. In fact, the 
critical velocity (Eq. 1) for the onset of the 
CVE is reduced by the factor y/(24.6/13.7)/4 ^ 
0.67 as compared to the pure hydrogen case. 
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However, the efficiency of the ionization is 
drastically reduced for two reasons. First, 
the resonant heating of the ionizing electrons 
depends on the growth rate of the modified 
two-stream instability which scales as the ion 
plasma frequency cUj. This frequency is much 
smaller for the heavier and much rarer he- 
lium component than for the hydrogen one. 
Consequently, the helium ions are expected 
to isotropize very fast and thus, cannot be 
ionized efficiently, i.e. e ~ 0.025, (see discus- 
sion by Raadu 1978; Formisano et al. 1982; 
Brenning 1985); Second, ionization is possi- 
ble only, if the Townsend criterion Uion > t~}^ 
if fulffiled where Vion and Te-i are the ioniza- 
tion frequency and the time an electron has 
to ionize an ion. Thus, the neutral gas den- 
sity has to exceed a critical value n„ > nc = 

[Te-iaion{Ee)^l2Ejrrie]~^ (where Ee is the av- 
erage energy of a single ionizing electron and 
aion is the ionization cross section). Whereas 
the ionization cross section for hydrogen and 
helium do not differ that much (see Mitchner 
& Kruger 1973) for our apphcations the par- 
ticle density n„ of the neutral helium is much 
lower than the one of HI. 

3. Ionization of the Cloudlets in the 
Magellanic Stream 

In the context of HVC's in the Galactic 
halo the CVE should work very efficiently. 
Concerning the high-velocity cloudlets in the 
MS, e.g. MS I- VI, several authors have estab- 
lished a set of values for the relevant plasma 
parameters as the particle density n, tem- 
perature T, and magnetic field strength B. 
Doppler measurements of Ha-emission lines 
show typical velocities v^^^ of the cloudlets 
in the range of 110 kms~^— 220 kms~^ (Weiner 
& Williams 1996; Murali 2000) for the cen- 
tral part of the MS and up to 350kms~^ for 



cloudlets at the tip of the Stream (Mirabel et 
al. 1979). The cloudlets show a rather low 
temperature of T^^^ 10^ K in a hot am- 
bient halo gas at a temperature of T^'^^° 
2 ■ lO^K (Weiner & Williams 1996). The 
particle density of the HVC's lies around 
^Hvc ^ 5 . 10-2 (Cohen 1982). How- 

ever, the particle density of the ambient halo 
gas cannot be determined so easily. Direct 
measurements of the halo density are not pos- 
sible so far since only column densities can be 
measured and the problem of the exact dis- 
tance of the MS still remains unsolved. The 
only way to get an estimate of the halo den- 
sity consists in the deduction of the density 
via the use of hydrodynamic models for the 
HVC's. From a rather simple drag model 
Weiner and Williams deduce a halo density 
of the order ^ 10"'' cm~^ while Murah 
finds a maximum value of 10"^ cm~^ from life- 
time considerations. The latter value, how- 
ever, does not seem very likely considering the 
fact that the particle density of the intergalac- 
tic medium already shows a minimum value 
of 10~^ cm~^. By assuming a pressure equilib- 
rium between the cloudlets and the ambient 
halo gas, Mirabel et al. find particle densities 
between 7 • 10~^ cm~^ and 7 ■ 10"'^ cm~^ de- 
pending on the distance of the MS. 
The above data are shown as single data 
points in Fig. 1 assuming a magnetic field 
strength in the halo of 3^°"^° 3yuG (Beuer- 
mann et al. 1985). For the mentioned halo 
plasma parameters the electron plasma beta /3e 
is much smaller than unity such that equa- 
tion (^ holds. The two dash-dotted lines 
in Fig. 1 mark the lower boundaries of the 
regions of efficient electron heating for the 
cases i^rci = 350kms~^ (upper line) and v^^i = 
220kms~^ (lower line). Obviously, the condi- 
tions for the onset of ionization of the HVC's 
boundaries caused by the CVE are satisfied 
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for all data points shown in the figure (see 
Eqs. 1 and 2). 

Consequently, ionospheres similar to the iono- 
spheric layers around planets and comets form 
around the HVC's. Due to these ionospheres 
more or less sharp boundaries form. Since the 
magnetic field is frozen in the halo plasma, 
it will be draped around the HVC's during 
their motion through the halo plasma. The 
convecting magnetic field generates currents 
in the conducting ionospheres that keep the 
field from penetrating through the clouds, at 
least initially. Depending on the conductivity 
and possibly local changes of the halo field 
direction some parts of the field may even- 
tually diffuse into the clouds. In any case a 
more or less perfect magnetic cavity is cre- 
ated. The draped field lines form a magnetic 
barrier. The magnetic field strength in such 
a magnetic barrier may well exceed the value 
of 10/iG as can be deduced from ram pres- 
sure equilibrium models (Kahn 1997). Closer 
to the centers of the HVC's ionopauses form 
where the internal ionospheric pressures bal- 
ance the incident pressure of the halo. Be- 
yond the magnetic barrier even a coUisionless 
bow shock may develop if the velocities of the 
HVC's exceed the local halo magnetosonic ve- 
locity. The described scenario for the forma- 
tion of boundary layers around HVC's is il- 
lustrated in Fig. 2. It should be noted that 
the whole scenario resembles very much the 
well studied interaction of the solar wind with 
unmagnetized planets (Luhmann 1986; Luh- 
mann et al. 1990) and comets (Formisano et 
al. 1982; McComas et al. 1987; Flammer et 
al. 1997). 

The Townsend condition for the produc- 
tion of secondary electrons that in course are 
heated and ionize the neutral gas further de- 
fines a maximum radius for the ionization re- 



gion, the so-called CVE-radius: 

i?cvE ~ l/n^'^'^aion ~ 10-' PC (5) 

where (Jion ~ 10~^^ cm^ is the ionization cross 
section (compare Flammer et al. 1997). 
The given CVEl-radius is very small com- 
pared to the the typical size of the cloudlets 
in the MS which is of the order of 1 kpc 
(Weiner & Williams 1996). This means that 
newly created ions at the cloud temperature 
should mainly be found in a boundary layer 
around the HVC. This however is consistent 
with observations of localized Ha-emission at 
the leading edges of the clouds (Weiner & 
Williams 1996). 

Furthermore, the cloudlets should be stable 
against full ionization by a Townsend dis- 
charge. 

Diffuse Ha-emission from the inner parts of 
HVC's is then due to remnants of the ioniza- 
tion front at earlier times. 

It is worth mentioning that the CVE also 
results in a braking of the cold HVC's. The 
ionized boundaries interact with the magnetic 
field of the halo that is mainly directed par- 
allel to the Galactic plane. The clouds are 
slowed down while the newly formed ions are 
accelerated. The overall effect consists in a 
slowing down of the ionization front. Ultra- 
violet absorption observations hint at a de- 
celeration of the HVC's as they approach the 
Galactic disk (Benjamin 1999; Danly 1989). 

4. Discussion 

In our paper we discussed the critical ve- 
locity effect as a possible explanation for the 
observed Ha emission at the leading edges of 
cold cloudlets in the MS. For suitable parame- 
ter regimes this effect can yield a rapid ioniza- 
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tion of neutral gas streaming through a mag- 
netized plasma. 

Our work has been motivated by intensive 
studies of the CVE with respect to observed 
ionization phenomena in laboratory and space 
plasmas. In laboratory plasmas ionization 
caused by the CVE has been clearly identi- 
fied in numerous experiments (e.g. Fahleson 
1961; Danielsson & Brenning 1975; Piel et 
al. 1980). The conditions, however, under 
which the CVE becomes sufficiently efficient 
in space plasmas are still under debate. 

Ionization due to the CVE has been clearly 
identified in context with a space experiment 
where an artificial cloud of barium was re- 
leased in the Earth's magnetosphere (Haeren- 
del 1982). Besides, ion production in the lo 
plasma torus around Jupiter and ionization 
in lo's ionosphere are interpreted as the con- 
sequences of the CVE (Cloutier et al. 1978; 
Galeev & Chabibrachmanov 1983). 

Further ionization phenomena in space which 
are discussed to be the consequence of the 
CVE are the space shuttle glow (Papadopou- 
los 1984) and, especially, ionization in cometary 
environments due to the interaction between 
neutral gas sublimented from the cometary 
nuclei and the magnetized solar wind plasma 
(Formisano et al. 1982). 

All these examples yielding evidence for 
the importance of the CVE refer to plasma 
- neutral gas interactions in laboratory and 
space plasmas which differ from a large scale 
astrophysical system hke the MS in several 
aspects. However, the presented first quan- 
titative studies of the CVE considering pa- 
rameter regimes being typical for cloudlets in 
the MS strongly indicate that the CVE can 
yield an efficient ionization of the neutral gas 
at the boundaries of cold cloudlets in the halo 
plasma. 



Thus, our kinematical studies have shown 
that the CVE has to be considered as a pos- 
sible explanation for the observed Ha emis- 
sions. To get deeper insight into the question 
of to what extent the CVE can contribute to 
the ionization of halo clouds detailed dynam- 
ical multifluid studies are required. 
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Fig. 1. — The CVE in dependence of the 
magnetic field strength and the particle den- 
sity (cf. Eq.dl) in double-logarithmic repre- 
sentation. The CVE can operate efficiently 
(e 1) for the parameter regime between 
the solid line and the dotted one (for the case 
Vrei = 10~^c) or the solid line and the dashed 
one (for the case Vrei = 10~^c), respectively. 
The data points represent the parameters for 
HVC's in the MS. 



Fig. 2. — The formation of boundary layers 
around HVC's. The CVE leads to the forma- 
tion of an ionospheric boundary. Probably a 
coUisionless bow shock forms as in the case of 
magnetized and unmagnetized planets. 
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